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The observation of the Kondo effect in meso- 
scopic systems under bias 1,2 has opened a new 
chapter in the physics of the Kondo phenomenon. 
Various types of dl/dV, where / and V denote cur- 
rent and source-drain (s-d) bias, respectively, line 
shapes have been measured by scanning tunnel- 
ing microscopy (STM) 1 ' 3-11 . However, explana- 
tion by single Fano line shape 1, is not relevant 
and even misleading. Here, we provide consistent 
explanations for various asymmetric and symmet- 
ric line shapes in terms of a microscopic theory 
that shows the creation of two resonant tunneling 
levels (RTLs) when bias is applied 17 . One side 
Kondo coupling between adatom and substrate 
does not create Kondo peak that appears only 
when the system has an overall Kondo coupling 
including both substrate and tip. The structure 
of an asymmetric line shape is mostly governed 
by the RTL peaks. Therefore, Kondo effect is 
negligible in most asymmetric line shapes. 

The dl/dV line shapes measured by STM 1,3 " 12 can be 
classified into four different groups. Typical examples of 
such groups are given as follows: (i) a Co adatom on a 
Au(lll) 1 or Cu(lll) 4 ~ 7 surface showing a dip near zero 
bias; (ii) a Co adatom on a Cu(100) surface showing a sig- 
nificantly asymmetric line shape 7,8 ; (iii) a Co adatom on 
an insulating CU2N film covering a Cu(100) surface show- 
ing a sharp Kondo peak with symmetric shoulders 9,10 ; 
and (iv) a system similar to (iii) with an Fe adatom 
showing a peculiar wedge and symmetric step-wise line 
shape 11 . 

Since the existence of two RTLs is the new discovery of 
this study and plays a crucial role in explaining the line 
shape, we discuss it first. The occurrence of two RTLs is 
a unique characteristic of nonequilibrium steady state of 
a mesoscopic system with two reservoirs. Because of the 
back and forth transits between the Kondo impurity and 
the left and right reservoirs, two more basis vectors are 
needed to describe the dynamics of the system under bias 
compared with the system of single reservoir. These two 
additional basis vectors yield two RTLs. In equilibrium, 
however, these additional basis vectors are not needed 
because the probabilities of back and forth transits are 
the same. Only relevant nonequilibrium theories that can 
properly handle the Kondo phenomena in a steady state 
can verify the existence of the two RTLs. This resonant 



tunneling is another type of coherent transport channel 
in addition to the Kondo channel. Interference between 
these two coherent scattering channels could yield Fano 
resonances 18 at the positions of the RTL because the 
RTLs are discrete states created within the continuum 
spectrum of the Kondo peak. However, only dispersed 
RTL peaks are observed in most STM line shapes in- 
stead of the sharp Fano interference pattern because no 
Kondo peak occurs except case (iii) mentioned above. 

We employ Tersoff-Hamann formula 19 that gives the 
differential conductance of an STM system as 



dl_ 

dV 



where T* denotes the featureless coupling function be- 
tween adatom and STM tip, and p„(V) is the local den- 
sity of states (LDOS) of the adatom at a steady state 
of bias V. The spectral function of an up-spin electron 
at the adatom is given by p"(u>) — (l/7r)Re[(M J T 1 ^ 
where 
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where to' = uj — e a — U(n a i) 18 . All the matrix ele- 
ments, except the 8 iJ-elements, have additional self- 
energy terms, /3 mn [iSg(cj) -I- iSpM] = 2/3 m „A, for a flat 
wide band, where A = (r s -fT')/4. We use A to indicate 
the unit of energy. 

We calculate the matrix elements of M r for the single- 
impurity Anderson model with two metallic reservoirs 
whose Hamiltonian is given by 

/c,er,a£ s,t 

+ Yj (^ka c aa c ka + ^ka c kl C aa)j 
/c,ct,ckG s,t 

where s, t, e a , and U denote the substrate, tip, energy 
level of the adatom, and Coulomb interaction, respec- 
tively. The detailed procedure to reach the matrix M r 
and its elements is given in Supplementary information 
of this paper. 
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The matrix M r contains all the information on the 
transport properties of the STM-Kondo systems. We dis- 
cuss the roles of 7 first and those of Uj± and the coeffi- 
cients j3 mn later. The upper-left 3x3 block describes the 
single-reservoir Anderson model between the substrate 
and the magnetic adatom, while the lower-right 3x3 
block represents another single-reservoir Anderson model 
between the tip and the adatom. The 2x2 blocks at the 
upper-right and lower-left corners couple the substrate 
and the tip through the adatom. Therefore, 7 ss (7tt) rep- 
resents the strength of the one-side Kondo coupling be- 
tween the substrate (tip) and the adatom, and 7 s t repre- 
sents the coupling mechanism between the two reservoirs. 
On the other hand, the antidiagonal element 77, where 
1j- = 7j+ = 7Ji gives the effect of steady state. 

The operators comprising 7, which are given in Sup- 
plementary information of this paper, give clearer inter- 
pretations for the role of each 7. We present pictorial 
descriptions of 7 in Fig. 1. Figure la shows the one-sided 
Kondo coupling of l ss r t t)> Fig. lb shows the difference 
between the leftward and rightward motions of a down- 
spin electron in 77; Fig. lc shows the overall Kondo cou- 
pling of 7 s t, i.e., an up-spin electron in both reservoirs 
has Kondo coupling with a down-spin electron passing 
through the adatom via exchanging or hopping together. 
Figures lb and lc also show the resonant tunneling of a 
down-spin electron if it does not make a Kondo coupling 
with an up-spin electron. 

In most of STM systems, no Kondo coupling exists 
between the adatom and tip, i.e., jtt = 0. In this case, 
the system does not have an overall Kondo coupling and 
1st = 7j- This implies that only resonant tunneling is 
allowed in Figs lb and c. However, when the system has 
the overall Kondo coupling, 7 st > 77 and both 7 SS and 
7tt are nonvanishing. In fact, the values of j s t an d 7j at 
a steady state cannot be determined by present theory. 
However, we find that choosing appropriate constants for 
7 s t and jj gives rise to the dl/dV line shapes that are 
well-agreed with experiments. 

Another important information is obtained from the 
matrix M r . The zeros of the determinant of M r in the 
atomic limit, i.e., A = 0, gives us the positions of five 
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FIG. 1: Motions of a down-spin electron that makes a Kondo 
coupling with an up-spin electron. The processes indicated by 
number two are those of exchanging and hopping together, 
a, 7 SS and ~/tt, b, jj, and c, j si . We do not show up-spin 
movements in b. b and c also show the resonant tunneling of 
a down-spin when Kondo coupling is not constructed. 




FIG. 2: Change in dl /dV line shape. The unit of the ordinate 
is (4e 2 //i)r*. We set [7 = 8. The unspecified parameters are 
the same as those used in the previous panel, a, /3u = 0.252, 
012 = 0.246, 022 = 0.256, and other mn = 0.25. j ss = 
0.12, 7 « = 0.1, 7 S < = 0.65, 7j = 0.35, Re[U a J+ ] = 2.4, other 
Re[Uj] = 2, and Im[U s /*] = 0. b, 0n = 0.255, i2 = 0.24, 
022 = 0.27, 7tt = 0, 7st = 7 j = 0.5, Re[U^] = 0.5, and 
Im[{7j + ] = 0.3. Dashed line: 7 SS = 0. c, ^ s t = 7j = 0.35 and 
lm[U S j+] = 0.8. 



peaks of the spectral function at the adatom. The five 
peaks are the Kondo peak at u' = 0, two Hubbard peaks 
at lj' ±U/2, and two RTL peaks at J = ±[{j 2 s + 
7 t 2 t )/2 + ( 7st - 7j ) 2 + Oil/- 2 )} 1 / 2 . The spectral weight 
of the Kondo peak can also be obtained from M r as 
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where the subscript 2R indicates a two-reservoir system. 
It is interesting to note that when 7 st = jj and j tt — 
0, which is the case of most STM-Kondo systems, Z 2 r 
vanishes and the gap between the two RTLs is V2^f ss . 
Since no Kondo peak exists in this case, the Kondo effect 
only appears as a separation of RTL peaks. One can 
see this case in the experiments on Co atom adsorbed 
on Au(lll) 1 or Cu(lll) 4 ~ 7 surface. On the other hand, 
when 7 s t > 77, i.e., ju 7^ 0; %2R does not vanish and the 
gap between two RTLs becomes wider. This implies that 
the two RTLs recede by level repulsion due to the Kondo 
peak at the Fermi level. One can see this case in the 
experiments on Co atom adsorbed on Cu(100) surface 
covering with Cu2N 9-10 . 

We now show our theoretical dl/dV line shapes ex- 
plicitly. Figure 2a, which is obtained by choosing 7 SS > 
7tt > 0, j st > 7j, and Im[C/j + ] = 0, clearly shows both 
the Kondo and Fano peaks. However, the Kondo and 
Fano peaks disappear and only the RTL peaks remains 
as shown in Fig. 2b if we choose 744 = 0, 7^4 = 7/, and 
Im[C/j + ] > 0, which are the usual case of STM-Kondo 
system. Figure 2c is obtained by using a larger asym- 
metry parameter Im[J7J + ]. The red lines in Figs. 2b and 
2c correspond to the dl/dV line shapes observed in a 
Co adatom on a Cu(lll) 4 and a Au(lll) 1 substrate, re- 
spectively. A small shift in the dip from zero bias stems 
from the nonvanishing asymmetry parameter. The con- 
ditions 7 s t = 77 and jtt = eliminate the Kondo peak 
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FIG. 3: Various dl jdV line shapes for U = 8. The unit of the 
ordinate and the presentation scheme are the same as Fig. 2. 
a, /3n = 0.265,^12 = 0.23, /3 22 = 0.28, and other /3 m „ = 
1/4. 7ss = jtt = 0, and lst = 7 j = 0.25. Re = 0.5, 

lm[Uj+] = 0.4, and other Im[i7 JT ] = 0. b, Im[U s J+ ] = 0.2. 
c, -y sa = 0.07 and Im[U s J+ ] = 0.3. d, /3u = 0.251, (3 12 = 
0.249 and (5 2 2 = 0.252. 7ss = 0.3, jtt = 0.1, 7s < = 0.5, and 
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7s t = U = 0.7, Re[U^] = 0.4, Re[U a J+ ] = 2, and Re[£/* + ] = 
1. Fermi level is shifted 0.1 to the right only in c. In h, 
two additional parameters have been used to take the spin 
excitation into account. 



in Fig. 2a according to equation |T]) and also changes 
the Fano resonance into a dispersed RTL peak because 
the Kondo spectrum, which is necessary for constructing 
Fano interference, has been removed. 

The blue lines in Figs. 2b and c indicate the line 
shapes without the Kondo coupling in both sides, i.e., 
7ss = 7*t = 0. The similarity between the red and blue 
lines indicates that the structure of the dl/dV line shape 
is mostly governed by the RTL spectrum. The Kondo 
effect between a Co adatom and a substrate appears as 
a slight widening of the dip structure and a shift in the 
minimum. The amount of widening may be considered 
as an energy scale of the STM-Kondo system with dip 
structure because no Kondo peak appears in this case. 
We measure the increase 6D in Fig. 2b, i.e., the differ- 
ence between the gaps in red and blue, as SD = 0.1D, 
where D is the full distance of the gap in red. Since D in 
rcf. 4 is approximately 20 meV, we obtain SD « 2.0 meV 
for the system used by Manoharan et al. 4 . This energy 
corresponds to a temperature of 23 K that is consider- 
ably higher than the temperature at which usual STM 
measurements are performed. In order measure the elon- 
gated gap, one must have a scheme to remove the Kondo 
coupling between the adatom and substrate. 

Highly asymmetric line shapes have been observed in a 
Co adatom on a Cu(100) surface 7-9 . We obtain similar 
line shapes observed by Neel et al. s in Fig. 3a and b 
by adopting 7 SS = "f t t = 0, which means that no Kondo 
effect is involved in obtaining these line shapes. However, 
a weak Kondo coupling between adatom and substrate is 



needed to obtain Fig. 3c, which is the line shape reported 
by Choi, et al. 9 for the same system. In conclusion, no 
substantial Kondo coherence exists in the system of a 
Co adatom on a Cu(100) surface 7-9 that shows highly 
asymmetric dl/dV line shapes. 

A remarkable change in the line shape occurs when 
the Cu(100) substrate is covered by an insulating cop- 
per nitride film. The line shape becomes symmetric, and 
a Kondo peak and side peaks 9 or steps 10 appear for a 
Co adatom. However, a peculiar wedge at zero bias and 
multiple steps appear for an Fe adatom 11 . Covering the 
Cu(100) surface with an insulating copper nitride film 
causes at least two effects. One is to strengthen the co- 
herence in axial direction connecting Cu in substrate, Cu 
in film, and Co adatom. The other is to suppress both 
fluctuation ((Sj^) 2 ) 1 ^ 2 and activity (j^f) on the sub- 
strate side. The former results in the Kondo coupling 
between the adatom and tip and the latter makes 7 SS 
and Re[(7^ + ] increase. The latter also causes a decrease 
m (Wal ) 2 ) 1 ^ 2 - As a result, all the imaginary parts of 
Uj' T vanish, and the line shape becomes symmetric. 

On the basis of this analysis, we obtain a symmetric 
line shape with the Kondo peak and two RTL peaks as 
shown in Fig. 3d. This line shape corresponds to the 
one reported by Choi et al. 9 . Selecting a relatively large 
Re[{7j+] changes the RTL peaks of Fig. 3d to shoulders, 
as shown in Fig. 3e, which has been reported by Otte 
et al. . An interesting change occurs when we simply 
choose 7 st = 7j and j tt = in Fig. 3e to remove the 
Kondo peak. A sharp wedge appears at zero bias as 
shown in Fig. 3f. The wedge reflects the trace of a miss- 
ing Kondo peak. Such a wedge has been observed by 
Hirjibehedin et al. by adopting Fe adatom instead of 
Co. We show in Fig. 3g that a step-like feature can be 
obtained by changing Re[L/J±]. An additional step ap- 
pears in the line shape of ref. 11. We understand that it 
stems from the spin-spin interaction at the position of the 
adatom. If we add the spin interaction at the adatom, 
the matrix M r becomes a 7 x 7 matrix and additional 
steps appear as shown in Fig. 3h. 

In conclusion, we have shown that two RTLs can ex- 
plain most of the experiments on STM-Kondo systems 
consistently. If the Kondo coherence exists only in be- 
tween a magnetic adatom and a substrate, the Kondo 
peak is not observed. The energy scale of the Kondo 
coherence in this case may be defined by the amount 
of gap elongation of the dip. Interestingly, the highly 
asymmetric line shapes are simply the RTL spectra with 
asymmetry. We showed that the Kondo peak and steps 
appearing in ref. 10 are attributed to the overall Kondo 
coupling and the suppression of the fluctuation on the 
substrate side. Finally, we understand that the sharp 
wedge for an Fe adatom reflects the trace of disappeared 
Kondo peak and the additional step stems from the spin 
interaction at the adsorbed atom. 
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Supplementary Information: 
The Procedure Reaching Matrix M r and the Expressions of /3„ 



The spectral function of a magnetic atom adsorbed on a substrate is given by p a }{u) = — (I / n)ImG^ a Juj) . We 

show below that the spectral function p a -\{u) or the retarded Green's function G^ a ^(u>) can be expressed by the matrix 

M" 1 given in the text of the main article. Since the retarded Green's function expressed in a resolvent form in the 
Heisenberg picture is written as 1 



aa\ V / 



where L is the Liouville operator defined by hO = [Ti.,0] for the Hamiltonian Tt, z = —iu) + + , c f indicates 
the annihilation of an up-spin electron at the adsorbed atom, and the matrix M is constructed by the elements 
M.ke = z8kt + (ee\iLek), where the vectors eg and et are the elements of a complete set of orthonormal basis vectors 
{ee\£ = 1, ■ ■ ■ , 00} spanning the Liouville space. The operator c a f is one of the members. The inner product is defined 
by (ek\ee) = ({efe,ej}), where the angular and curly brackets denote the expectation value and anticommutator, 
respectively. 

Therefore, the first step to obtain the retarded Green's function G^ a ^(ui) is to find a complete set of orthonormal 
basis vectors of the Liouville space. In ref. 2 given below, we have shown a basis spanning the reduced Liouville space 
that is effective for describing the dynamics of the system in the Kondo regime as 

(c| T , Sn al c s kv Sj+ L s c ah Sj-fcai, c ah 6j-fc ah 6j+{c ah 6n al c{ v c* feT ) 
with normalization factors ((^n a j) 2 ) -1 / 2 and ((Sj^f ,t ) 2 )^ 1 ^ 2 for the corresponding basis vectors, where SO = O — (0), 



■T s ,t 



3ai'~ = W Sfe(^feo c ai c fc| T ^fea' c k{ c a|): where (i) is deleted for j a ^ , and k = 1, 2, • • • , 00 indicates the quantum 
states of the substrate or the tip. 

We construct the full matrix M using these basis vectors in an order given above. Then, M is given by the following 
nine-block matrix, 



M ss M as 
M = I M sa M a M ta 
M at M tt 



(SI) 



where the blocks M a , M as and M at , and M sa and M ta are 5 x 5, 5 x 00, and 00 x 5 matrices, respectively. Since 
no direct coupling exists between the substrate and tip, null blocks occur at the two corners. The block M ss ( tt ) is 

constructed by the basis vectors eff and 5n a ic^\ It is written as 



M 



ss(tt) — 



M11 
Mn 



where 



/ z + iex •••0 \ 



Mn = 



•• 

V •••0 z + ie x ) 
On the other hand, M sa = — M.\ s and the 5 x 00 block M as has the form 



M„ 















c s ._ a .+ c ., c* ._ 



where the columns C/c a , C*'._ , and C*'.+ have elements iVj* a -, — (£ s '*Vfc a )*, and — (£<J~ 5 ' Vfc )*, respectively, where 



= (l/2)[(iK i5 jJ^](l - 2n aT )) + i(l - 2(n al ))Ul''%({6j: i '' t r)({5n al r)] 



-1/2 
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The first term is purely real and the second term is purely imaginary. The detailed process of obtaining ^J s ' t T4 a from 
— (iLcfcf Sn a j.|c T<y^), which is the element of M sa , is given in Appendix A below. 
The 5x5 block M a at center is given by 
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where Cj = uj — e a — U(n a ]). 7 and Uj^ are given by 



lss(tt) 



ist = C£i(v k * a c s k1 +K aC ^ T )ci T b7/,^])K(^-/) 2 >((^ a + i t ) 2 )r 1 



/2 



7J T = (J2^ V ^h +^a4T)ci T b a T /,^]>[((^) 2 )((^a T l t ) 2 >]" 1/2 , 



and 



U 



iU_ 
T 



(KjJ^jl - 2n a] )) + - 2(n al )) 

((Ci s '*) 2 ) 1/2 



The detailed process of calculation of the matrix elements is given in Appendix A in which the Hamiltonian for the 
single reservoir Anderson model is used. 

The infinite-dimensional matrix M of equation (SI) can be reduced to an equivalent 5 x 5-dimensional matrix using 
Lowdin's partitioning technique 3,4 . For this purpose, we consider an eigenvalue equation MC = 0, where C and 
are column vectors. We partition the column vector C into three parts, i.e., C = (C s C a C*) T , where T means the 
transpose and C s , C , and C t correspond to rows (c| T , <5n aJ c| T ), {6j+*Cai, Sj^c^, c a t, 8j~*c a i, 5j+*c a i), 
and (Sn a ic k p c^), respectively. After eliminating C s and C t , we reach an equation (M a — MjjM^Mdj — 
M ta M M 1 M at )C a ee M r C a = 0. Then, we have the expression of the reduced matrix M r given in the text, i.e., 



/ — iuj + iT,ii(u>) 
7 ss +iE 2 i(w) 

-UT- 

-l st + ffi 4 i(w) 

V -l.i- + »S 5 i(w) 
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— TP* 
U J+ 
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where iT. mn {uj) = /3 m „(i£g(u) + = 2/3 mn A for a flat wide-band, where A ee (A s + A*)/2 = (T s + T*)/4. 

A is used as a unit of energy. The coefficients (3 mn are symmetric. They are given by fin = |£~ s | 2 , /3 22 = |£+ s | 2 , 

/?33 - 1, /?44 = |CT, /?55 - |CT, &4 = /?12 - &4 = "C*. 016 - C^C, 025 = C^C, and 

/?45 = £j~**£ a The coefficients /3 TO „ are complex when m ^ n, while /3 mm are real. The detailed expressions of /3 m „ 
are given in Appendix B. 



Appendix A 

Calculation of the matrix elements of M for the single-impurity Anderson model with one reservoir is shown below. 
Since calculation for the elements of M ss ( tt ) is very simple, we skip it. 



(1) Matrix elements of the block M sa : 

Nontrivial elements of M sa arc — {{iL{ck-\5n a {), 6j^, c' a A), i.e., 
{{i[H,Ckt6n al ], i^4 T }) = ({i[H,Ckt]5n al +c feT i[fl',n J,^4 T }) 



The first term 
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(— iek)({ckfSn a i, Sj^c' a ^}) = (— iefe)([^oj.)i^Jcfe|C^|) should vanish because the up-spin dynamics of this form 

is not allowed. The third term also vanishes because {{c^j~y, Sj^c^ }) = (Jai^Jaii^h c ai}) = ^- ^hc socon d term, 
however, must be calculated rigorously, since it is a hybridization term that is related to the Kondo process. The 
operator {c a ^Sn al , Sj^c^} is expanded as c a] -[Sn al , Sj^c^} + [Sj^c^, c a ^]8n ai + 2c a] Sj^ lC ^6n al . First two terms 

are rewritten as c^[5n al , Sj^c^j = c„ T [n al , ]4 T = c oT 4 T [n oi , jj] = (1 - n oT )[n ai , j^J and 

[^.4t> c <>t] ( * , VI = ^'04 [ c if c aTl^"-ai- = ^ijC 1 - 2c aT c l T )^ai, while the third term is cancelled by the last 
expression -2c oT ^4 T 5n oi . Hence, -({iL(c feT 5n oi ), j^4 T }) = iV fea ((l - n oT )[n ai , (fy'^J + 5^(5n ai ) = 
*V fca ((i - n aT ) [n oi ,^J + ±Ki, J a T J + ^5n oi ) = ^((1 - 2n oT )[n oi ,^J + {5n oi ,^}) = ^{((1 - 
2na T )[nai,^]> + (l-2<n oi ))0^>}. 

Therefore, the matrix elements — ~^~~~~~~~~7^~77~~~^'~ = Vfc £f , where 

l|c fcT <5n ai ||x||c aT <5j^|| sos< " 

£ = §{<t(l - 2n aT )[n ai ,iJ]) + i(l - 2(n ai ))(.^)}((fc ai ) 2 )-V2 ( ( (5j T )2 > -i/2. 

(2) Matrix elements of the block M a : 
(2-1) Diagonal elements: 

By using the commutator expression i[H, c a |] = — *Ek^ka c kT — * e aC a | — iJ7c a fn a |, the diagonal elements of the 
block M a are given by 

(A) - ({i [H, c a] ] , 4 T }) - i E k ^ka ({ c kT . 4t » + lt * ({ C «T . 4 T » + * U{{c a] n ai , 4 T }) 
= ie a + ill ({c aT n a ;, 4 T }) = + ^ <{c a t, c it}"-al) = + iU ( n <4>> and 

(B) -({i[^,c aT ^],(c aT ^)t}) = -({i^c^]^,^^^}) - ({c aT i[^,^],^4 T }). The first term of (B) is 
rewritten as 

({{[H^a^Sj^Sj^}) = ({HEkC c kt - ieaCai - iUc a] n al )8j^ (c aT <y^)t» 

= -^{{c^Sj^Sj^c^}) -iU {{c^na^j^Sj^c^}) . Applying the decoupling approximation, n ai <5.^ = (n a {)5j* v 
to the [/-term above gives rise to a form of squared norm for the first term of (B), i.e., ({i[H, c a ^]Sj^, Sj^c^ }) = 
{—ie a — iU (n a i)] x | |c a f<5j~J | 2 . The second term of (B), however, cannot be written in the form of a squared norm, 
because [H, j^jj cx j^. Therefore, we neglect it. 

(2-2) Matrix elements Ujt ■ 

These arc given by the inner products -({iLc ah Sj^c 1 ^ }) , i.e., -{{i[H, c |], Sj^c^}) 
= ie a (Sjl)({c ah 4 T }> + iU{{ca^n al ,Sj^}) + 1 E k V&Qj*) 4 T }) 

= if/({c oT n oi ,^4T})- Howcvcr : ({caT^^^lt}) = ( C «T [^M, #^4?]) 

+ {[Sj^i 4t> c a|]^ n oi) + (2 c aT < ^Ja F |4T^ ri ' a i)' wnere the operator in the first term is rewritten as 

CailSna^Sj^c^] = c a i[n al , jJ]4 T = c oT4tKi'iai] = C 1 ~ "el) Ki , jj] , while the second term is given by 
([ 5 jai c lv Ca ^ Sria ^ = (^[^T'^Tl^al) = ( S 3ai 5ria ~ ( 2S jal c ^ c U Sn »l) ■ The last one cancels the third term 
of ({catSna^Sj^c^}) above. Hence, ({c^Sria^Sj^c^}) 

= ((1 - n a] )[n al ,j^ l } + Sj^Sn^) = ((± - n„ T ) [n al , jJJ + ±[<5n ai , (5jJJ + <5j'^(5n ai ) 
= i((l - 2n oT )[n oi ,iJ] + {5n oi ,<yjJ}) = ±{<(1 - 2n oT )[r M ,j£]> + (1 - 2(n ai ))(.^ )}. 

Therefore, the matrix elements f7j T that are defined by — ^^g^^i/V^ are gi yen by 

Ujt = ^((1 - 2n oT )[n oi ,jJ]) + i(l - 2<n ai ))0 a T )]<(^)V 1/2 

If we define Rc[C/ JT ] = {/ /2r, we get r = 2 by the atomic limit analysis for the two-reservoir Anderson model 5 . Then, 
we have a relation, 1/r = — 2n a f )[n |, j^)((5j^) 2 ) _1 / 2 ; which will be used in Appendix B. 

(2-3) Matrix elements J ss (tt) '■ 

The inner products — ({iL(cof5j^), 5j^4t }) gi yes r i se to the expressions for 7 ss (tt)- Since 
({^(ca^j^),^^^}) = ({i^CaiJJj^^j^c^}) + ({ca^lHJ^^Sj^c^}), the first term can be expanded 
as ({«[-ff,c aT ]^,(5j^ci T }) = -^({ca-iSj^Sj^c^}} - iUdc^na^j^Sj^c^}} - i E k ^ko({ c fcT^. 6 jtl c U})- 
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The first and second terms of this expression should vanish because they are not written in the form of a squared 
norm. Since the third term describes the Kondo process, we treat it rigorously and have the following expression: 

-*£k^<{cAV44 T }> = -zEk^a^aWit - <y£^ fcT 4 T }) ' 

The second term of ({iL(c a ^Sj^), Sj^c^}} , i.e., {{c a -[i[H, j^J, Sj^c^}), is rather complicate because it contains 
the commutator [H,j^]. This commutator is expanded as i[H,j^] = ±e a j^ T («) J2k e fe(^4a4j. c ol ± V ka c ll c u) ± 
U n atjai ^~ W i^VkaVtacl^cei ± V fe * a V^* 4j. c fci ) , where (i) is applied only to the lower signs and the prime in sum 
denotes k ^ I. The first two are cancelled if we assume e a ~ and the [/-term can be neglected in the Kondo regime. 
The last term describes the round trip of a down-spin electron between the adatom and the substrate (tip), which 
is equivalent to (? a j) 2 - Since we consider only single trip of a down-spin electron in this work, we neglect the term 

l • [rr -Tl rrn , ■ , , A ,21/12 . . , -{{iL(c a1 5j^ ),5j± c+ }) i V V k * Q <c fcT c+ [jj j± ]) 

involving LH, 7 T , . Inus, tne matrix eiement M a is given by — , 1 . l , 1 = , ,' , t = ±7 

for the single-reservoir Anderson model. For the two-reservoir Anderson model, similar calculation gives rise to 

7as(«) = 



((Jj -«(t) ) 2 > l/2 <(5j +-(*) ) 2 > l/2 

Appendix B: /3 mn 

We obtain the expressions of (3 m n in terms of the definitions and the expressions of £J given in Appendix A. 
The coefficients [3 mn in front of the self-energy function, i.e., iY, mn (u)) = /3 TO „(iSo(w) + H2q(uj)), are symmetric in 
exchanging their indices. We have obtained an expression 1/r = — 2n a ^)[n a \,, ij])((^ij) 2 } -1 ^ 2 m Appendix A. 
By using the operator identity [n j.,J^] = TU^, onc can obtain a different expression for r, i.e., 



1/r = [±0j s, *> T 20-J s,t n„ T >]/V<WaT t ) 2 >. 



which will be used in the expressions of Re[/3 m „]. The final forms of Re[/3 m „] are given by using {(Sn a i) 2 ) w 1/4 and 
t = 2. Since 0'7^) < and > 0, an inequality, Re[/3i 2 ] < Re[/3n] < Re[/3 22 ], is applied to the substrate 

side in which Kondo coupling exists. As for Im[/3 m „], we neglect all of them because Im[/3 mn ] contain the difference 
of relative fluctuations, which we assume they are the same. 

We obtain the expressions of (3 mn in terms of given in Appendix A as follows: 
Real parts of /3 mn : 

Re[/3n] = ReKrSl 



1 , (i-2(»^))^o a -;)^ 



+ 



<(-5i.7) 2 > 



(0„Y> - 2 ^Tia + /)) 2 + (! " 2(n ai )) 2 ( Ja 7) 2 

_^ [1 (i-2(» ai » 2 (j-°; 

4<(5«ai) 2 ) ~ 4 + <(«"/) 2 > 



[4((fa^) 2 )]~ 



Only the index s changes to t for Re[/?5s]. 

Re[/3 22 ] = Refc+'tf] = [«J„7> - 2<n„ T j 7» 2 + (1 - 2<r M » 2 0#> 2 ' 

(1-2(»^)) 2 Q+;) 2 ] J _ [ t | (l-2(n tti )) 2 (j+ s ) 



««7) 2 > 



[4((fa,j) 2 )]- 

<w7) 2 > 



+ 



4((<5n ta ) 2 ) 



<(«7) 2 > 



Here, again, the index s changes to t for Re[/?44]- 
Re^ 12 ] = Re[C*C] = MP21} 

(0 a 7> - 2<n oT j+'))(-0 o 7> + 2<«aTJ„7» +(! - 2<^i)) 2 a 7)0'7) 

(i-2(» ai » 2 (j7)<j+;) 



! (l-2<n ai » 2 0-/><j+ s ) 
77 H ■ 



4<(<5"ai) 2 ) 



v / «^7) 2 >v /((,5 ^ S)2) 
Re[/3 14 ] = Rc[£7& + ] = Re[/3 4 i] 

(0„7> - 2<n oT i7))(-0 o 7) + 2<n-TJa7» +( X ~ 2 M) 2 a 7>0 a 7> 



! (l-2<n„ i » 2 07>07> 



4<(5"ai) 2 ) 



1 (l-2(n^)) 2 Q a 7)Q-2-/) 



-\/<(^7) 2 >\/<(«7) 2 > 



Re[/3 15 ] = Re[C^r] = Re[Asi] 

= "((,+/) -2<n oT j+«))(0-+j t > -2<n oT j+*)) +(1 - 2<n oi >) 2 0'7>0„7> 



y/(w:n a )y/(w:;) a ) 



(l-2(r Iai )) 2 (j- s )(j- t ) 



y/^j-p')y/((Sj-*)') 
Re[/3 24 ] = Re[f+*& + ] = Re[/3 4 



4((Sn a j)*> 



i (i-2(» ai )) 2 a 7)0 a -;) 

4 V^T^V^O^ 



(-0„7> + 2<n oT j o 7))(-0 o - J t > + 2(n aT i-/}) +(1 - 2(n ai )) 2 (,+/)(,+*) 



[4((fe ai ) 2 >]- 



! (l-2<n ai » 2 0-+ s >(j+*> 
_| 



4<(<5n ai ) 2 ) 



1 (l-2{»^)) 2 0- + 3 )Q-+ t ) 

4 



VWr^vTTTT 2 } 

Re[/? 25 ] = Rc[e+*C] - Re[/3 52 ] 

= "(-07/) + 2<»aTJ o "i S >)(0'#> - 2 <n aT j+ T >) +(1 - 2(n ai )) 2 a + s )0 a -/} 



[4((fo al ) 2 )]- 



1 (l-2<n„;» 2 <j+;)<J- t ) 

-r2 + ~ ' 



4((<5n ta ) 2 ) 



1 (l-2(»^)) 2 Q a + 3 )0 a 7) 

4 v««jj;) 2 >v«0 2 > 



Re[/? 45 ] = Rc[e t + *C] = Re[/3 54 ] 

- '(-07/) + 2(n aT i-/))(0 a + /) - 2<r M j+'» +(1 - 2(n ai )) 2 (. ? +/)( Ja 7 



v7(*77V^77 2 ) 



^((fa^) 2 )]- 1 



\Z<Wai) 2 >\/<wr/) 2 > 



1 , (l-2<n a; ))-<jffl<j-') 

^2 T 



4<(5« ai ) 2 ) 



(l-2(n ai )) 
4<(<5»o i ) 2 ) 

(l-2(n ni )) 1 
4<(<5«ai) 2 ) r 



v 7(<o77 2 )a/7 , o77 2 ) 

arts of /3 mn : 
?7*£7 

Q-+;>-2(n aT3 +;) o+p 



Imaginary parts of (3. 

Im[/3 12 ] = ImlC*^] 



077 077 



1 , (l-2(n^)) 2 ( J +/)0 a - t ) 

4 y<^^^> 



-o a ;)+2<"a T j a ;) o a ;) 

V<(«£) 2 > y/(W:-)') 



v7(<o77 2 ) v /((^7/) 2 > 

Im[/3 14 ] = Im[C*e t + ] 

Q-+;)-2(n aT3 +;> 



_ (l-2(» al )) 
~ 4((5n oi )-) 

_ (l-2(n ai )) 1 
" 4((5n a j)*> r 



-07/)+2(n aT j- f ) p--;> 



V((^7;) 2 > V ((5 77 2 ) 



v7(<o77 2 ) V((*jrj*) 2 > 

Im[/3i 5 ] = Im[£-*£ t ~] 

0+;)-2(n aT j+;> p--*) 



- ( 1 -2{"ai)) 

_ 4((5n ai )*> 
_ (l-2(n oi }) 1 



Jm\R 1 - (i-2("<a)) 
lm^ 25 j - 4((faaJ 2 ) 



v7w:7 2 ) v'<wr J t ) 2 > 

077 _ 07/) 



vW7 2 > v^u^:;) 2 ) 



Q-+ t )-2<n aT j+ t ) Q-;) 



(l-2(n ai )) i 
4<(<5n ai ) 2 ) r 



V<( 5 v) 2 ) v^c^)J' 

07j. s )+ 2 ("aTJ7j. s ) 077 

07/) 077 

a/((«7-/) 2 ) A /««7+;) 2 ) 



Im[/3. 



15 J 



(l-2(n aJ .)) 



(1-2(m)) 1 

4<(<5"ai) 2 > T 



077 



-077+ 2 <"»t47) 077 
v/« 5 77 2 > y/(W;{)') 



077-2(n aTJ 77 077 
v^wtT 2 ) vTTovTTh 



07)- 2 <"°t77 Qa + 7 



(77 



>/((^7/) 2 ) V ((5 77 2 >. 



and 



Im^ 24 j - 4((feai) 2) 



-o a 7>+2(» aT j a 7) (j+i) 



-07)+2<»a T 377 07 s ) 



(l-2(n al » 1 0+7 
4<(<5"ai) 2 > r 



V<Waj") 2 > \/<( 5 77 2 ) 
07) 



V« 5 77 2 > V« 5 77 2 ) 



V<( 5 7/) 2 ) V« 5 77 2 ) 
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